Effect of a magnetic field on the two-phonon Raman scattering in graphene 
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We have studied, both experimentally and theoretically, the change of the so-called 2D band of 
the Raman scattering spectrum of graphene (the two-phonon peak near 2700 cm"'^) in an external 
magnetic field applied perpendicular to the graphene crystal plane at liquid helium temperature. A 
shift to lower frequency and broadening of this band is observed as the magnetic field is increased 
from to 33 T. At fields up to 5-10 T the changes are quadratic in the field while they become linear 
at higher magnetic fields. This effect is explained by the curving of the quasiclassical trajectories 
of the photo-excited electrons and holes in the magnetic field, which enables us (i) to extract the 
electron inelastic scattering rate, and (ii) to conclude that electronic scattering accounts for about 
half of the measured width of the 2D peak. 

PACS numbers: 63.22.-|-m, 63.20.Kr, 78.30.Na, 78.67.-n 



I. INTRODUCTION 

Graphene, a two-dimensional plane of carbon atoms 
arranged in a honeycomb lattice, has stimulated many 
experimental and theoretical efforts to study and un- 
derstand its exotic electronic and optical properties^. 
Graphene can be produced by mechanical exfoliation of 
bulk graphite^, by the decomposition of a SiC substrate 
at high temperatures^ or by epitaxial growth on differ- 
ent metal surfaces^. Graphene on SiC is well adapted 
for optical measurements as samples are of macroscopic 
size and lie on a substrate which can be insulating, and 
hence, transparent in a wide, from far-infrared to visible 
spectral range. 

Raman scattering is a key element to graphene studies 
as it is a reliable and non-destructive technique to estab- 
lish the monolayer character of a graphene specimen^^S or 
the negligeable interlayer coupling in multilayer epitaxial 
graphene (MEG) samples^. This relies on the analysis 
of the shape of the observed 2D band feature (the two- 
phonon peak around 2700 cm~^, also referred as D* or 
G') which is the second-order overtone of the D band. 
While defects (in a general meaning - impurities, inter- 
face roughness, etc.) are needed for the D band to be 
observed, the 2D band is always seen and it even appears 
as the most pronounced feature in the Raman scattering 
spectrum of graphene. Because of the doubly-resonant 
nature of the D band*'^ and the fully-resonant nature of 
the 2D band,^*''^^ they both show a dispersive behavior 
with the excitation energy, what allowed to tracing the 
phonon band structure of different carbon allotrope o^^i^"^ . 
In this paper, we show that a strong magnetic field ap- 
plied perpendicular to the plane of a graphene crystal 
significantly affects the 2D band feature. Increasing the 



magnetic field up to 33 T, we observe a simultaneous red 
shift of the 2D band energy and a strong broadening of 
the peak. 

In the following, we show that the evolution of the 2D 
band in magnetic field can be understood and modelled 
quantitatively as a pure orbital effect of the magnetic field 
on the intermediate states of the Raman process, i. e., 
electrons and holes. The field induces circular orbits to 
the photo-generated electron-hole pairs, and this modifies 
the momenta of the optical phonons emitted during the 
Raman scattering process. This leads to a red shift and to 
a broadening of the 2D band feature as the magnetic field 
is increased, an effect that we observe experimentally and 
describe theoretically below. We note that this effect is 
not specific to a graphene monolayer as the two essential 
ingredients are the fully-resonant nature of the 2D peak, 
and the free motion of electrons in the plane. Thus, it 
should also be characteristic for the multilayer graphene 
and graphite. 



II. SAMPLES AND EXPERIMENT 

Raman scattering measurements were performed in the 
backscattering geometry at liquid helium temperature 
and in magnetic fields up to 33 T applied perpendicu- 
lar to the plane of the graphene crystal. A Ti:Saphire 
laser tuned at an accurately controlled wavelength of 
720 nm (corresponding to the incident photon energy 
hujin — 1.722 eV) was used for excitation. Optical fibers, 
both with a core diameter of 200 /im, were used for exci- 
tation and for collection. The resulting laser spot on the 
sample had a diameter of ~600 /im with a typical power 
of -100 mW. 



The samples presented in this study are 
by the thermal decomposition of the carb 
SiC substrate^. Despite the fact that these 
tain ^70 layers, simple Dirac-like electronic 
as the ones found in exfoliated graphene 
persist in these highly graphitized structure 
the particular rotational stacking exhibited 
cent graphitic planes^. Dirac Fermions ii 
structures have been evidenced by magneto — 
by magneto-transmissioi*i^ and, more recer 
ning tunnelling spectroscopy^^ and by ma^ 
scattering experimentsiS.. Apart from the 
terface layers which are highly doped due t 
transfer from the SiC substrate^^, the ma 
which we probe with Raman scattering, are ( 
with densities as low as 5 x 10^ cm~^ and 
high as 250000 cm2(V • s)"\ as deduced frt 
transmission experiments^^. MEG sample 
cially the highly graphitized specimens, are 
nous on the scale of our laser spot and c 
Bernal-stacked inclusions as observed thr 
Raman scattering experimentSi^. 

Typical Raman scattering spectra of the 2 
sured at B = 0, 10, 20 and 30 T are present □ 
At zero field, the 2D band is observed as a s' 
metric feature composed of two Lorentzian C' 
The main contribution to this feature is 
2646 cm~^ with a full width at half maximi 
of 26 cm~^ and the second, weaker, contrib 
tered at 2675 cioT^ with a FWHM of 31 cm" 
ously observed on similar samples^, the ene 
the 2D band is observed is increased, for i 
tation wavelength, with respect to the enei 
the 2D band is expected for exfoliated gra 
The origin of this effect, probably due to 
type of stacking of the graphitic planes and 
electronic interaction between these planes, i„ ^ ..^^„ 
ter of debatc^'^'^^. The origin of the double component 
2D band feature observed in our experiment, different 
than the one observed in bulk graphite, is still not clear. 
It could be due to Bernal-stacked inclusions under the 
large laser spot. We note that similar line shapes have 
also been observed in micro-Raman scattering measure- 
ments on epitaxial graphene samples on SiC and were 
interpreted as revealing graphene domains with different 
amount of strain32.. When a magnetic field is applied per- 
pendicular to the graphene crystal plane, both of these 
two components show a red shift and a broadening. We 
present in Fig. [5] the evolution of the Raman shift and 
of the FWHM of both components of the 2D band with 
magnetic field. From i? = to 33 T, a shift of ^ 8 cm^^ 
and an increase of 20% of the FWHM is observed. 



III. DISCUSSION 

From the very beginning, we would like to stress that 
the observed effect is quite unlikely to be related to the 
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FIG. 1: (Color online) Experimental Raman scattering spec- 
tra measured at T = 4.2 K in the 2D band range of energy 
at different values of the magnetic field (black solid line), two 
Lorentzian fits (red solid line) and independent components 
of the Lorentzian fits (black dotted lines). The vertical black 
dashed line is a guide for the eyes. 



modification of the spectrum of the emitted phonons 
(e. g., by electron-phonon interaction). Indeed, elec- 
trochemically top gated graphene structures allowing to 
tune the Fermi level up to 800 meV (carrier density up to 
^ 4x10^'^ cm~^) have been produced recently.^- Ra- 
man scattering experiments in these highly doped struc- 
tures showed that the evolution of the 2D band central 
frequency with doping is moderate and monotonous, and 
can be described by a change of the equilibrium lat- 
tice parameter. Modification of the phonon dispersion 
due to the coupling of phonons with low energy elec- 
tronic excitations across the Dirac point is noticeable 
near the K point of the phonon band structure (Kohn 
anomaly), while the phonons involved in the 2D band 
are too far from the K point to be affected by the Kohn 
anomaly. This argument remains valid even in magnetic 
fields used in our experiment. Indeed, the magnetic field 
mixes wave vectors on the scale of the inverse magnetic 
length, 1/Ib = y/eB/{hc) w 0.022 A" ^ at B = 33 T. The 
wave vector of the phonons responsible for the 2D peak. 
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FIG. 2: (Color online) a) FWHM and b) Raman shift of the 
low-energy component of the 2D band as a function of the 
magnetic field, c) FWHM and d) Raman shift of the high- 
energy component of the 2D band as a function of the mag- 
netic field. Solid red lines are calculated with our model (see 
text). 



measured from the K point, is q ~ 0.24 A at 1.7 eV 
excitation. Thus, g 3> 1/Ib and the Kohn anomaly re- 
mains inaccessible for these phonons even at our highest 
fields. This is in striking contrast with the main first- 
order Raman feature, the so-called G band due to the 
doubly degenerate i?2g optical phonons at the T point of 
the phonon band structure. These phonons are directly 
in the region of the Kohn anomaly, and their frequency 
and lifetime are noticeably affected even by a moderate 
gate voltage i^^i^^ In strong magnetic fields, the G band 
exhibits a pronounced magneto-phonon effect, as the 
phonon state is modified by the coupling to electronic 
inter-Landau-level transitions i^iS^ 

In the following, we show that our observations are due 
to the finite curvature of the trajectories of the photo- 
excited electrons and holes in the magnetic field. The 
main idea is illustrated in Fig. [3] In the semi-classical 
real-space picture of the Raman process,^- the incident 
photon creates an electron and the hole with opposite 
momenta at an arbitrary location within the laser spot. 
They subsequently propagate along the classical trajec- 
tories, and emit phonons. If they meet at some other 
location, again with opposite momenta, they can recom- 



bine radiatively producing a scattered photon. In the ab- 
sence of the magnetic field, the trajectories are straight 
lines, so that in order to meet at the same point with 
opposite momenta and contribute to the Raman scatter- 
ing signal, the electron and the hole must necessarily be 
scattered backwards during the phonon emissioui^ This 
fixes the phonon momentum hq (measured from the K 
or K' point) as q = p + p' , where hp = hu!in/{2v) and 
hp' = h{uji-a — 2a;ph)/(2w) are the electronic momenta 
(also measured from the Dirac points) before and af- 
ter the phonon emission, cjjn is the excitation frequency, 
Wph is the phonon frequency, and v is the electronic ve- 
locity (the slope of the Dirac cones). A slight spread 
in the phonon momentum due to the quantum uncer- 
tainty, \q — p — p'\ ~ 27/1; (27 being the electron inelas- 
tic scattering rate), gives a contribution to the width of 
the 2D peak, ^ 2^{v-p\Jv), where fph = duj-p\Jdq is the 
phonon group velocity^S 

In a magnetic field, the electron and hole trajectories 
are no longer straight lines but, because of the Lorentz 
force that acts on charged particles in a magnetic field, 
they correspond to circular cyclotron orbits. As a result, 
one can see from Fig. [3] that (i) phonons with smaller 
momenta, q = pcos(p + p'cosip', can be emitted, and 
(ii) since each phonon can be emitted at an arbitrary in- 
stant in time, the length of the arc describing the electron 
trajectory is random [not exceeding the electron mean 
free path v/ (27)], and so are the angles (p, ip' . Note that 
the frequencies of the two emitted phonons remain equal 
even with an applied magnetic field due to the symmetry 
under C2 rotation around the axis perpendicular to the 
crystal plane. Tuning the magnetic field is, in this sense, 
equivalent to changing the resonant conditions of the Ra- 
man scattering process at a fixed excitation wavelength 
and allows exploring part of the phonon band structure 
closer to the K point. Fact (i) results in the overall red 
shift of the Raman peak, while fact (ii) introduces an 
additional spread in q, and contributes significantly to 
the broadening of the peak as observed through Raman 
scattering measurements. 

To produce the observed effects, the electron and the 
hole do not have to complete a full orbit. The probability 
of this latter event vanishes as g-^'^f^i^f/^) ^ where R — 
p{eB / hc)~^ is the cyclotron radius. This implies that the 
separation between the electronic Landau levels, hujc = 
hv/R, may still be smaller than their broadening 2/17 
and, in our experiment, this is still the case even at B = 
30 T. In other words, even magnetic fields far from the 
quantization limit can give a noticeable effect. A closely 
related effect has been discussed in the context of density- 
density response of a degenerate electron gas in a non- 
quantizing magnetic field^. 

To obtain some more quantitative information about 
the behavior of the 2D peak in magnetic fields, we cal- 
culate the Raman matrix element A4{q) for a given 
phonon wave vector g in a clean graphene monolayer. 
It is convenient to perform the calculation in the coor- 
dinate representation, analogously to the calculation for 
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FIG. 3: (Color online) Schematic of the electron and hole 
motion during the Raman scattering process. The lightening 
represents the incident photon which creates the electron-hole 
pair. The solid arcs denote the propagation of the electron 
and the hole in the magnetic field. The fiash represents the 
radiative recombination of the electron-hole pair. The dashed 
arrows denote the emitted phonons. 

the D peak in the vicinity of an edge^ (details of this 
calculation are presented in the Appendix). An essen- 
tial ingredient of the calculation is the semiclassical elec- 
tronic Green's function for graphene in a magnetic field 
which was calculated in Ref. l32l Under the assumption 
LOc ^ (7,'^ph) ^ (^in, the result of the calculation can be 
represented as: 

M{q) cx J rf^7Je-W«-2p)+27H.-.b/(i2fl^)].^ ^ 


= 2(i7r)3/2 A A. \Ai\u) -iAi{u) Bi(u)] , (1) 

For vanishing magnetic fields, the cyclotron radius R — > 
00 and the relation A4{q) oc {q — 2p — li'^ jvf'l'^ , which 
is Eq. (64) of Ref. [2^, is recovered. For finite values of 
the magnetic field, the matrix element is expressed in 
terms of Airy functions Ai(M) and Bi(u)^. Note also 
the similarity with the expression for the polarization 
operator in Ref. HO. 

Since the density of the final phonon states is practi- 
cally (7-independen1*2i, \M.(yi\^l v + ^ j (2v^\^)\^ describes 
the Raman scattering intensity as a function of $7, the 
Raman shift measured with respect to the center of the 
2D peak at zero field. To illustrate the change of the 
peak shape in magnetic fields, we present in Fig. Ilja) 
|7W(wi„/'y r2/2t;ph)P for B = 0,10, SOT, considering 
Swin = 1.7 eV, hv ^ 1 ■ k [v = 1.06 x 10^ cm/s), 
'L'ph/i' = 50 cm~^/eV, hrf = 27 meV. This value for 
7 is chosen in order to reproduce the observed shift of 
the peak maximum for increasing magnetic fields. How- 
ever, the resulting FWHM for the 2D band is about 
twice smaller than the one observed experimentally at 
zero field. Most likely, this indicates the presence of an 
additional broadening mechanism besides the electronic 
scattering. As we cannot determine the nature of this 



mechanism, we model it by introducing phenomenologi- 
cally an additional broadening of the 2D band through 
a convolution of with a Gaussian curve of width a: 

hom^ / -^e M{ — + - . 3 

— CO 

There is no simple relation bet ween the FWHM of the 
original \M\'^ peak {8j{vph/v)V2'^/^ — 1 at zero field), 
the FWHM of the broadening Gaussian (2o-\/ln4), and 
the FWHM of the resulting peak. The result of the con- 
volution for a = 9.5 cm~^ is shown in Fig. 

In the following, we assume 7 and a to be indepen- 
dent of the magnetic field. This approximation is rea- 
sonable as long as the magnetic field is non- quant izing.— 
At B = 30 T, electrons with energy e — 0.85 eV mea- 
sured from the Dirac point, have a cyclotron frequency 
ujc = {v'^h/e){eB/hc) of about 27 meV, so the strongest 
magnetic fields in this experiment seem to be close to 
the limits of validity of this approximation. Beyond this 
regime. Landau quantization of the electronic spectrum 
should manifest itself as oscillations of the peak intensity 
due to periodic modifications of the resonance conditions 
with increasing magnetic field, while in the experiment 
no such oscillations are observed. 

The results of this calculation in terms of Raman shift 
and of FWHM of the peak, as described by Eq. (O, 
are presented as red solid lines in Fig. [2] (a)-(d) for the 
two components of the observed 2D band feature. The 
following parameters were adjusted: (i) the central fre- 
quencies of the two components at zero field, giving the 
overall vertical offset for the curves in Fig. [2] (b), (d); 
(ii) a — 9.5 cm^^ and a — 11.9 cm~^ were taken for the 
low-frequency and high-frequency components, respec- 
tively, in order to reproduce the FWHM at zero field, 
in combination with (iii) /17 = 27 meV which determines 
all four slopes in Fig. [5] (a)-(d). Clearly, without precise 
knowledge of the origin of the two 2D band components 
we cannot account for the slight difference in their zero- 
field widths. Nevertheless, the fact that the calculation 
reproduces the red shift of the 2D band energy as a func- 
tion of the magnetic field which is quadratic for fields up 
to 5-10 T and linear at higher fields, with a single value 
of 7 for both components, shows that the electrons have 
the same dynamics in the parts of the sample, responsible 
for the two components. 

Besides, the deduced value of hj ~ 27 meV is in good 
agreement with that deduced from the doping depen- 
dence of the 2D peak intensity in exfoliated graphene^i^ 
(given the energy dependence of 7 and the fact that the 
latter measurements were performed at higher excitation 
energy, the electron scattering in those samples is some- 
what weaker than in ours). We also note that the found 
value = 27 meV is in reasonable agreement with the 
line width of electronic transitions in high magnetic fields 
measured in this range of energy on similar samples 
This fact is remarkable because those measurements were 
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FIG. 5: Integrated intensity of the 2D band low energy com- 
ponent (black dots) and of the 2D band high energy compo- 
nent (open circles) as a function of the magnetic field. 

performed in a quantizing magnetic field, and a priori the 
scattering rate does not have to be the samei^ 

From Eq. ([IJ it can be seen that the integral 
J \M{q)\^dq does not depend on the magnetic field. (It 
is sufficient to integrate over q first, and then over z; 



the magnetic field enters only through R, which drops 
out.) This means that under the assumption of a con- 
stant phonon density of states, the frequency-integrated 
intensity (the area under the peak) of the 2D band should 
not depend on the magnetic field. The experimental in- 
tensities for the two components are plotted on Fig. [S] 
Only the integrated intensity of the high-frequency com- 
ponent is field independent, while an overall decrease of 
^ 35% of the integrated intensity of the low-frequency 
component is observed over the range of B between 
and 33 T. Does it mean that the phonon density of states 
decreases stronger for the lower- frequency component? 
Again, without precise knowledge of the exact nature of 
the two components it is hard to give an explanation for 
their different behavior. 



IV. CONCLUSIONS 

In conclusion, we have studied, both experimentally 
and theoretically, the evolution of the 2D band of MEG 
structures in intense magnetic fields. We observe a red 
shift and a broadening of the 2D band as the magnetic 
field is increased. We have modelled this effect using a 
semi-classical picture in which the Lorentz force induced 
by the magnetic field on the photo-created electron-hole 
pairs, curves the carriers trajectories. This leads to a 
decrease of the emitted optical phonon momentum and 
a broadening of the 2D band as observed through Ra- 
man scattering spectroscopy. This model enables us to 
extract the value of the electronic scattering rate. From 
this we conclude that about half of the 2D band width 
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at zero field is due to electronic scattering while the ori- 
gin of the other half remains to be clarified. Finally, we 
note that the observed effect of the magnetic field on the 
2D peak is not specific to monolayer graphene; it should 
be analogous for multilayer graphene and graphite. 



where indices i,j = x,y,z label the Cartesian compo- 
nents, and the momentum is defined as 



e r 

V Irl 



■A(r). 



(A4) 
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Appendix A: Quasiclassical calculation of the 
two-phonon Raman matrix element in a uniform 
magnetic field 

1. Classical action for a Dirac particle 

The two-phonon Raman matrix element is given by a 
loop of four electronic Green's functions. The latter will 
be taken in the quasiclassical approximation. The most 
important ingredient of the quasiclassical Green's func- 
tion is the classical action. Thus, we first discuss the 
classical action, then we give the explicit expression for 
the electronic Green's function, and finally, we perform 
the calculation of the Raman matrix element. In fact, the 
first two steps have been already made by Carmier and 
Ullmo2^, but we include them here for the sake of com- 
pleteness of the presentation, and to fix the notations. 

We start with the classical equation of motion of a 
charge in the magnetic field B: 



dp e 

— = - V X B 

dt c 



(Al) 



For a Dirac electron with energy e the velocity v(t) = 
vn{t) and momentum p(f) = (e/i;)n(t) are expressed 
in terms of a unit vector, |n| = 1. [If instead of elec- 
trons with positive and negative energies one prefers 
to work with positive-energy electrons and holes, then 
p(i) = |e/w|n(t), so the action given below should be 
multiplied by sgn e] . Eq. ()Al|) can be obtained by varia- 
tion of the action functional 



S[v{t)] 



j dt |r| + - A(r) • f 



(A2) 



keeping the ends ri,r2 of the trajectory fixed. Indeed, 
upon integration by parts 



5S 



dt 



e d Xi e 
V dt |r| c 
P2 • 5t2 - Pi • 5yi, 



dxi dxj 



Sxi 
(A3) 



If we now define the function S{r,r') as the action on the 
classical trajectory, corresponding to the motion from r' 
to r according to Eq. (|A1I) . we obtain VS* = p at the 
end of the trajectory. Thus, this function satisfies the 
Hamilton-Jacobi equation 



V5(r,r') - - A(r) 

c 



-V'5(r,r')- -A(r') 
c 



(A5) 



Also, dS/de is equal to the time it takes to go along the 
trajectory. 

In a uniform magnetic field, n is precessing with a con- 
stant frequency ui = — (ei?/c)(v^/e), and the trajectories 
are circles of the radius R = v/\u}\. Two particular points 
r and r' can be connected either by two trajectories cor- 
responding to short and long arcs (plus an integer num- 
ber of full rotations), or by no trajectories at all if the 
distance between the points is greater than the circle di- 
ameter, |r — r'l > 2R. Let us assume B to be along the 
z axis and choose the gauge A(r) = {B/2)[ez x r]. The 
action along the short/long arc is given by 



5±(r,r') = ^[^±(r,r')+sin^±(r,r')] 
2v 



eB_ 
"2^ 



[r X 



d+ (r. r') = 2 arcsin ■ 



2R 



??_ (r, r') = 27r — 2 arcsin 



|r — r 
2R 



(A6) 
(A7) 
(A8) 



Here i?±(r, r') is the angular size of the short/long arc. 
We also introduce i?2^(r, r') = 27rj+i9±(r, r'). We denote 
by n-t the unit tangent vector to the short/long arc at 
the point r (Fig. [5]): 



n±(r,r') 



|r — r 
r — r' 



sm sgn uj 

2 ^ 

^?±(r,r') 



(A9) 



The tangent determines the direction of the kinematic 
momentum: 

VS'±(r,r') - -A(r) = -n±(r,r'), (AlO) 

c V 

-V'5±(r, r') - -A(r') = n^(r, r'). (AH) 

c V 
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R 



FIG. 6: (Color online) Two classical trajectories connecting 
the points r, r'. 



2. Quasiclassical Green's function 

Now we discuss the Green's function for the Dirac 
equation, 



e + ih'jf^ sgne + v S- (ihV + -A) \ G(r,r' 



= ]l(5(r-r'). (A12) 

Here S = (T,x,'^y) are the pseudospin matrices in the 
sublattice space, 1 is the unit matrix, and 7^ > is 
the half of the electron or hole inelastic scattering rate, 
introduced phenomenologically (see Sec. IVC of Ref. [23 
for the discussion of its role in Raman scattering). The 
factor sgne corresponds to the chronologically ordered 
Green's function. 

The quasiclassical Green's function is represented as a 
sum over all trajectories of the classical motion from r' 
to r, i. e., sum over the short and long arcs, "±", and 
over the number of rotations j — 0,1, . . .: 



G(r,r') 



E 



GL''(r,r') 



(A13) 



G2'^(r,r') = e±"/4 



sgne. 



,\/ihvy 



X i^n 



±(r,r' 



n^(r,r') 



g-i(eB/2nc)[rxr'], ^ 



X e 

^ g-)?2''(r,r')7.-R/i'^ 



(A14) 



The imaginary part of the argument of the exponential is 
just the classical action discussed in the previous subsec- 
tion. Replacement e — |e| is required by the analytical 
properties of the chronologically ordered Green's function 
(otherwise it would correspond to the retarded Green's 
function). The damping factor e~'''^^^'" is determined 
by the total length of the corresponding arc -dR, includ- 
ing 2TTjR from j full rotations. The eigenvector ip^ of 
(n • S)i/'n = V'n sgne, |n| = 1, can be written in terms of 



the polar angle (pn of the unit vector n (we use the basis 
where S^;, Ey are represented by Pauli matrices). 



■011 



-i¥>„/2 



V2 V e*'^"/2 



sgne 



(A15) 



Thus defined, V'n is not a single-valued function: when 
n is rotated by 27r, ipn acquires a minus sign. We can fix 
the sign by requiring that at r — >■ r', tpn+ — — '0-n_ = 
'i/'r-r'i V^n- — 4'-n+ — V'r'-ri and they evolve continu- 
ously during the motion along the circle. 

We make several further remarks concerning the above 
expression for G(r,r'), Eqs. (|A13|) . (|A14|) . 

(i) For r — >■ r' G(r, r';e) reduces to the expression for 
the quasiclassical Green's function for B = 0, 



Go(r-r';e) = - 



j7r/4 



1 sgn e 



r — r 
|r — r' 



Ji\t\/h~',,)\r~r'\/v 



27r|r-r' 



(A16) 

valid at |e||r - r'\/{hv) > 1. 

(ii) The difference between the exact Green's function 
and G(r, r';e) given by Eqs. (|A13p . (|A14|) is of higher 
order in h. This can be checked explicitly by calculating 
the gradients: 

[n^ X dr], [n^ x dr']. 



V sin i?j 



Rsm-d± Rsm'd± 
1 cost9± r — r' 
R cos(i9±/2) |r-r'|' 



1 



-Is ^ 



(S • V)^^^ = - 



(n_ • I])V'n_ 



(A18) 

(A19) 
(A20) 



2i?sini9+ ' 
so that at r 7^ r' 

[e-|-i7£Sgne-t-iftt;S • V]G^^ (r, r'; e) = 

[n+xS],Gfi!")(r,r';e). (A21) 



ihv 



i?cos2(i?+/2) 

It is proportional to the eigenvector corresponding to the 
energy — e; thus, it can be compensated by a correction 
to the pre-exponential factor which does not vanish when 
acted upon by e -I- ivH ■ (VS* — ieA/c), and thus is of the 
next order in h: 



SG^i\r,r';e) = 



ihv 



2ei?cos2(i9+/2) 



[n+xS],G(^'(r,r';e 



(A22) 
aram- 

eter hv/{eR). However, it has a stronger divergence at 



(5G^' is smaller than G^'' by the dimensionless param- 
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|r— r'l — >■ 2R; this is the manifestation of the usual break- 
down of the quasiclassical approximation in the vicinity 
of the classical turning point. 

(iii) To fix the relative phases of C?^' and G^'' we 
choose a given circular trajectory and consider the evolu- 
tion of G(r, r'; e) as r moves along the circle. Namely, as 
r passes the turning point at |r — r'| = 2R, G'f\r,r';e) 

should transform into g'^I^v, r'; e). In the vicinity of the 
turning point, |r — r'| ~ 2R{1 + z), z ^ 1, we have 
'&± « tttV^- Thus, i?± can be viewed as two branches 
of the same analytical function, and — >■ ■)?_ when 
z makes a circle around z = in the complex plane. 
Since ■d± + sm'd± w tt =F (— 2z)^/^/6, in order to get a 
decaying exponential in the classically forbidden region 
|r — r'l > 2R (the positive semiaxis of z), the circle must 
be counterclockwise. Thus, 



-in/4 



sin-i? 



sint? 



0)i 



(A23) 



precisely as in Eq. (|A14|) . 

(iv) To match the phases of G^"* and c'f^^), we con- 
sider r — > r' and require that G^^ (r, r'; e)+G^^^^ (r, r'; e) 
satisfies the homogeneous Dirac equation. Indeed, the 
(5- function term is produced by G^^"^ component, see re- 
mark (i). At r r' the Green's function satisfying the 
homogeneous Dirac equation can be can be constructed 
in the i? — >■ limit. Let us introduce Go(r — r'; e), the 
Green's function with analytical properties opposite to 
those of Go- 



Go(r~r';e) = 



t/4 



-1 sgne 



27r|r- 



(A24) 



It also satisfies Eq. (jA12[) . as seen from the antisymme- 
try of the Dirac operator with respect to the simultaneous 



change e 



-J> -e, I -I, 7e -5> 



-7e. Thus, Go(r-r';e) 



Go(r— r'; e) represents the sought solution, corresponding 
to the flux of particles focusing at the point r'. Hence, 
the correspondence is g''1\y r';e) ~ — Go(r — r';e) 

(converging wave), G^^^' (r r'; e) ^ Go(r — r'; e) (out- 
going wave), which produces precisely the combination 
of signs as in Eq. (|A14[) . 

Remarks (iii) and (iv) can be restated in very simple 
terms. Upon every half-rotation sin d changes sign, so the 
square root produces a factor of e"/^ . In addition, upon 
a full rotation acquires the Berry phase e*'^. Hence, 
the sum over j in Eqs. (jA13l) . (jA14[) reduces to a simple 
geometric progression, 



j=0 



1 



on{i\e\/h-2'r)R/v ' 



(A25) 



which determines the poles corresponding to ireR/hv be- 
ing an integer multiple of 2tt, i. e., Bohr-Sommerfeld 
quantization rule. It gives the exact expression for the 
Landau levels, Icn/vl = ^y2n\ehB /c\. 



3. Two-phonon Raman matrix element 

Here we consider the matrix element (q; lyJm, V'out) 
for the transition from the initial state, corresponding 
to the incident photon with frequency uj-m polarized at 
the angle ip-m to the x axis, and no phonons, to the fi- 
nal state, corresponding to the scatterd photon with fre- 
quency Win — 2a;q polarized at the angle (^out to the x axis, 
and two phonons with momenta q, — q and frequencies 
cjq (q is measured from the K or K' point, and for the 



phonons in the two valleys the frequencies w, 

due to the time-reversal symmetry; we denote w, 
and omit the valley index) . Then, for unpolarized excita- 
tion and detection, the frequency-resolved intensity /(f2) 
is given by 



dfin dipout f d^q 



27r 27r J (27r)2 



X |X(q; <fin, 'Pout)\^S{n ~ 2wq). (A26) 



The matrix element is given by the loop of four Green's 
function with electron-photon and electron-phonon ver- 
tices (here we write it in the coordinate representation)^: 

7W(q;(/Jin,<Pout) oc J ^ J (fr2 (fvi (fro (fr[ 

X Tr {(Sj, cos ipout + Sy sin tpout) 
X G(r2,ri;e + Win/2 -cjq) 
X S,e-^i"-iG(ri,ro;e + Win/2) 
X (Sj. cos (fin + T^y sin ipin) 
X G(ro,r;;e- Win/2) S^e'i"-! 
X G(ri, r2;e- Win/2 -fwq)}. 

(A27) 

We will assume 7tJi„/2-R/i' ^ 1 (non-quantizing field), 
then the probability of a half- or full rotation is exponen- 
tially small, and only the G^'' contribution to the Green's 
functions remains. The integral over e is dominated by 
small |e| ~ 7wi„/27 so the energy arguments of the Green's 
functions can be assumed to have definite signs. Thus, 
the imaginary part of the exponent can be written as 

S'(r2,ri; Win/2 - Wq + e) - qri + 5'(ri, ro; Win/2 -f e) 
-hS'(ro,r;;win/2-e)-|-qr; 

+ 5(r'i,r2; Win/2 -wq + e), (A28) 



where the action S is that given by Eq. (IA6p . with the 
subscript omitted and the energy argument explic- 
itly introduced. 
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The spatial integration is performed in the station- 
ary point approximation. Namely, in the whole 8- 
dimensional space (fq, ri, r'^, we separate a manifold 
on which expression (jA28l) is stationary. As VS* is just 
the classical momentum, this manifold corresponds to 
joining the classical arcs connecting the pairs of points 
in order to satisfy momentum conservation at each point 
To, ri, r']^, r2. Then the integration over the deviations 
from this manifold is performed in the Gaussian approx- 
imation, while the integration over the manifold itself 
has to be done more carefully. This procedure is fully 
analogous to that used in Sec. VIB of Rcf. 31 for the 
edge-assisted Raman scattering. 

The integration over e will be performed by expanding 
the actions to the linear order, e. g. 



5'(ri,ro; Win/2 + e) 
S'(ri,ro; Win/2) - 



9S'(ri,ro; u;i„/2 + e) 



de 



(A29) 



e=0 



and neglecting the e-dependence of the pre-exponential 
factors. Then the integration over e gives a J-function, 
S(tQi+ti2 — t2i' —tiio), where Uj is the time of travel from 
the point i to the point j according to the classical equa- 
tions of motion. This (5- function simply expresses the fact 
that the electron and the hole have to travel for the same 
amount of time before the radiative recombination. 

Let us assume the phonon momentum q to be along 
the y axis. Out of eight spatial integration variables two 
correspond to translations of the trajectory as a whole. 
They contribute to normalization, but as we are not inter- 
ested here in the overall prefactor in J\4{q), they can be 
discarded. In the remaining 6-dimensional space we in- 
troduce three coordinates yi, ip, ip' which parametrize the 
stationary manifold, and three deviations 5xi,5R,5R' . 
The integration variables are parametrized as (see Fig. [7]): 



To 



+ {5R cos (/3, 5R sin Lp) , 
(0,yi) + (5x1,0), 
(0,-j/i) + (-fei,0), 



(A30) 
(A31) 
(A32) 

-R' cos ip' - ^J{R')^-yf, R! sin ^'^ 

+ {-6R'cosp',6R'&mLp'), (A33) 

where R = p/{eB/c), R' = p'/{eB/c), p = ujin/2v, p' = 
(win — Wq)/2w. If we iutroducc ipi — arcsin(yi /i?) , p'l = 
arcsin(?/i/i?'), then momentum conservation reads as 

psinpi = p' smp[, pcospi +p' cosip'i — q . (A34) 



This ensures that the expansion of Eq. (jA28[) in the de- 
viations does not have linear terms. 

Now we have to expand Eq. (jA28P to the second order 
in Sxi, SR, 6R'. If we denote 



r r 



(A35) 



P H 




FIG. 7: (Color online) Classical trajectories determining the 
stationary manifold. 



the second derivatives of the action can be written as 
dS dS 



dxidxi 



dx[dx'^ 



eB 
2^ 
dS 



—SiSj tan — t + (Sij — SiSj) cot 



,(A36) 



eB 



dxidx'^ 2c 



{5ij - SiSj) cot 



(A37) 



where Cxv 



1; ^XXZ 



0. As a result, 



for the quadratic part of the action AS* we have 
eB sin 2ipi {SR ~ Sxi cot ipi sin ip)'^ 



AS 



2c sin((/3i — ip) sm{ipi + ip) 
eB sin 2<p[ {6R' — 6x'i cot (p^ sin ip')'^ 
' 2^r 



sm{(p[ — ip') sin(v5i -|- (p') 



eB 



— (cot (pi + cot (p[){Sxi) 



(A38) 



After integration over dR,6R',6xi, performed in the 
Gaussian integration, we are left with three spatial in- 
tegration variables (fjip'^yi, and the energy variable. As 
mentioned above, the latter gives the temporal 5-funtion, 

S{toi + ti2 - t2v - ti'o) = 1 5{Rip + i?V') , (A39) 

which also lifts the integration over ip' . Let us take 
the limit Wq <C Win, and neglect the difference between 
R and R' and between ipi and ip'i. Choosing (pi as the 
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integration variable instead of yi , we obtain 

t/2 (^1 



X 



M(q) oc / Kcosipidipi / di/? 

J J sin2(^i 



^ 22pi?(v3i+sin cpi COS cpi) — 2zQi?sin c^i— 47i?cpi/v ^ 

Ac A 

X COS^ (ySi C0s((^in — </3) cos((/?out — f) ■ (A40) 

The first integration is simply R cos tpi dtpi = dj/i , the 
factor -^/i? tan sin 2(^1 comes from the gaussian inte- 
gration, and the cosines give the angular dependence of 
the electron-photon [cos(<pin,out— <(5)] and electron-phonon 
(costal) matrix elements. Integration over tp gives 

7r/2 



M{q) (X / i?^/^ difi \J sin Lpi cos^ Lpi x 



X e 



2zpi?(c^i +sin ipi cos ip\)—2iqB, sin ipi —A'yRipi/v 



COs(yin " lygout) 

(sinv?i)/<^i 



-I- COS^Si COs((pi„ -I- (/Sout) 



(A41) 



short arcs, so that Lpi and |2p — gr| ^ p. Thus, we ex- 
pand the exponent to ip\, while in the prefactor we keep 
the leading term at ifi — 0. In the experiment described 
in the paper, the excitation and detection are unpolar- 
izcd, so we simply omit the the expression in the square 
brackets, which describes the polarization dependence. 
Denoting R(pi — z, we arrive at 



M{q)o, / (A42) 



which is Eq. (1) of the main text. The integral is calcu- 
lated using the relations^^i^ 



dt 

2^3/2 ] 7f 



oo 

1 f dt 
2^ J Vt 





cos xt 



sin xt 
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= Ai2(x), (A43) 
= Ai(a;) Bi(a;XA44) 



Taking advantage of the limit 'jR/v 3> 1, we can focus on 
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